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ABSTRACT 



o 

■ We calculate the ion-ion correlation effect on the neutrino-nucleus scattering 
in supernova cores, which is an important opacity source for the neutrinos and 

^ ' plays a vital role in the supernova explosion. In order to calculate the ion-ion 

■ correlation effect we use the results of the improved hypernetted-chain method 
calculations of the classical one-component plasma. As in the preceding studies 
on this effect, we find a dramatic decrease of the effective neutrino-nucleus scat- 
tering cross section for relatively low energy neutrinos with E,j < 20 MeV. As a 
matter of fact, our calculation shows a much more dramatic reduction of the ef- 
fective neutrino-nucleus scattering cross section for the low energy neutrinos with 
Ey < 10 MeV than the results of Horowitz. Therefore, the ion-ion correlation 
effect will be more important than has hitherto been recognized. We present an 
accurate analytic fitting formula that summarizes our numerical results. This fit- 
ting formula will facilitate the application of the present results to the supernova 
explosion simulations. 
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1. INTRODUCTION 

In a recent important paper, Liebendorfer et al. (2003) have made thorough comparisons 
of various neutrino radiation-hydrodynamical codes of the supernova explosion caused by the 
delayed neutrino-driven mechanism. They have carefully assessed the magnitude of the ion- 
ion correlation effect on the neutrino-nucleus scattering in the supernova cores and have 
found that this effect causes an appreciable difference in the supernova core collapse results. 

Historically, Itoh (1975) pointed out the importance of the ion-ion correlation effect on 
the neutrino-nucleus scattering in the supernova cores for the first time shortly after Preed- 
man's (1974) theoretical discovery of the coherent aspect of the neutrino-nucleus scattering 
and its importance to the supernova problem. Itoh (1975) used Hansen's (1973) rather hm- 
ited Monte Carlo results on the classical one-component plasma then available and showed 
the appearance of the dramatic reduction of the effective neutrino-nucleus scattering cross 
section caused by the ion-ion correlation effect for relatively low energy neutrinos. Bowers & 
Wilson (1982) subsequently published an analytic fitting formula for the ion-ion correlation 
effect that is based on the work of Bond (1980). Bruenn (1993) discussed the consequence of 
the ion-ion correlation effect on the neutrino-nucleus scattering in supernova cores in great 
detail. 

More recently Horowitz (1997) presented a detailed fitting formula for the ion-ion cor- 
relation effect based on his Monte Carlo calculations. Bruenn & Mezzacappa (1997) and 
Liebendorfer et al. (2003) have made a detailed assessment of the importance of the ion-ion 
correlation effect on the neutrino-nucleus scattering in the supernova explosion using the 
results of Horowitz (1997). Since this effect plays an important role in the supernova explo- 
sion, we wish to carry out an independent calculation of the ion-ion correlation effect on the 
neutrino-nucleus scattering and make a detailed comparison with Horowitz's results in this 
paper. 

The present paper is organized as follows. Physical conditions relevant to the present 
calculation are made explicit in § 2. Calculation of the ion- ion correlation effect on the 
neutrino- nucleus scattering in supernova cores is presented in § 3. An analytic fitting formula 
that summarizes the numerical results is presented in § 4. Concluding remarks are given in 
§5. 

2. PHYSICAL CONDITIONS 

We assume that a supernova core consists of atomic nuclei {Z, A) and degenerate elec- 
trons. At the densities and temperatures relevant to the supernova core, these nuclei interact 



-3- 



with each other strongly by the Coulomb interaction. The conventional parameter that char- 
acterizes such strongly-interacting classical ions is (Itoh et al. 2003) 

where aj is the ion-sphere radius 

1 /3 

ai = 0.7346 X 10-^2 (^^) cm, (2) 

Z being the atomic number of the nucleus, A being the mass number of the nucleus, Tio 
being the temperature in units of 10^°K, and pi2 being the mass density in units of 10^^ g 
cm~^. 

In this paper we restrict ourselves to the case in which electrons are strongly degenerate. 
This condition is expressed as 

1/2 



T < Tp = 5.930 X 10 



9 



[l + 1.018{Z /Af/'pf] -1 



K, (3) 



where Tp is the electron Fermi temperature and pe is the mass density in units of 10^ g 
cm~^. We impose a further condition that the ions are in the hquid state (Nagara, Nagata, 
& Nakamura 1987), 

r < 172 . (4) 
The conventional r^-parameter for the electron hquid is given by (Itoh et al. 2002) 

rs = = 1.388 X 10-^ f-^y^' , (5) 

-nal rie = 1 , (6) 
where rie is the electron number density. The electron Fermi wave number is expressed as 

/Z \ 

kp = 2.613 X 10^M-pi2j cm-^ (7) 

We have the relationship 

kpaj = [yT' Z^'- (8) 

When electrons are extremely relativistic, the Fermi-Thomas screening wave number is 
expressed as (Flowers & Itoh 1976) 



kr- ' 



vpT 



kp V TT 137.036 



0.0964 . (9) 
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Therefore, for the extremely relativistic electrons, one obtains 

kpTtti = 0.1850Z^/^ (10) 

In Figure 1 we show the density-temperature diagram of the pure ^^Fe plasma that illustrates 
the physical conditions described in this section. Here one should note that we assume the 
ion-ion correlation effect is caused only by the Coulomb interaction between the nuclei. At 
extremely high densities p > 10^^ g cm~^ one should also take into account the nuclear 
(strong) interaction to calculate the ion-ion correlation effect accurately. 



3. ION-ION CORRELATION EFFECT 

Bruenn & Mezzacappa (1997) have given detailed expressions for the ion- ion correlation 
effect on the neutrino-nucleus scattering in the supernova core. We refer the reader to this 
very comprehensive paper for the detailed discussion of this effect. 

It is well-known in liquid metal physics that one has to include the liquid structure factor 
to obtain the effective cross section for the X-ray scattering in the liquid metals (March 1968). 
This condition equally applies to the neutrino-nucleus scattering in supernova cores, as the 
nuclei are strongly correlated with each other by the Coulomb interaction. The strength of 
the ion-ion correlation is measured by the Coulomb coupling parameter F in equation (1). 

The effective transport cross section for the neutrino-nucleus scattering is given by 
(Freedman 1974; Itoh 1975) 

^) = al^A^El{l + cose) {1- cose) S{k), (11) 

where e is the neutrino laboratory scattering angle, Oq is a numerical parameter of the 
theory of the weak interaction, G is the Fermi coupling constant, A is the mass number of 
the nucleus, and S{k) is the liquid structure factor. The momentum transfer k is related to 
the neutrino energy E^, by 

k — — — sm-. (12) 
he 2 ^ ^ 

We refer the reader to Bruenn & Mezzacappa (1997) for more detailed expressions for the 
cross section. 

For the uncorrelated ions with S{k) — 1, we have the angle- integrated cross section 

< >uncorr — 3^0^^^ ■ i^^) 
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<S> 



Therefore, we can express the ion-ion correlation effect S{k) ^ 1 by writing 

4 

< cr >corr = -al—A^El <S>, (14) 
^ J dcos9 (1 + cos^) (1 - cos^) S{k) 

= 12 /" x^{l-x'^)S{x)dx, (15) 
Jo 

_ ajk ajk , , 

^ = 2E,ai/hc ^ ' ^ ' 

_ E^ai _ E^ai , , 

^ " ~ 197.3 MeVfm ' ^ ^ 

So far we have considered only one nuclear species {Z, A). However, in the real supernova 
explosion more than one nuclear species exist in supernova cores (Liebendorfer et al. 2003). 
Here we give a prescription to calculate the neutrino mean free path i due to the neutrino- 
nucleus scattering including protons: 



0-1 



= J2nj<aj>, (18) 



where Uj is the number density of ions of j-th species and < aj > is the neutrino-nucleus 
scattering cross section by j-th species ions (nuclei). In calculating < aj > for the neutrino- 
nucleus scattering by j'-th nuclear species {Zj,Aj), one uses the mass number Aj and the 
angle-averaged ion- ion correlation factor < 5* > in equation (14) corresponding to the F- value 
(Itoh et al. 1979) 

r, = = 0.2275 ^ Pi2 y ^ , (19) 

where Xj is the mass fraction of the i-th nuclear species (Zj, Aj) and Og is the electron-sphere 
radius 



3 \ 1/3 / Q \ 1/3 



A-nrieJ \47r Hi Zi 
Ue being the electron number density. 



(20) 



For the ^^Fe plasma with Z—26, we have kpTai—O-^'^^ from equation (10). Hubbard 
& Slattery (1971) have carried out a Monte Carlo simulation of the ions screened with the 
dielectric function of the electron liquid. They have shown that for kpTCLi < 0.5 the results 
that take into account the electron dielectric function are close to the results corresponding to 
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the ions imbedded in the negative background. Therefore, we are justified in using the hquid 
structure factor corresponding to the ions imbedded in the negative background (classical 
one-component plasma). 

In this paper we use the liquid structure factor of the classical one-component plasma 
calculated by Ichimaru, lyetomi, & Tanaka (1987) using the improved hypernetted-chain 
scheme. They have presented a numerical tabic for r=2, 5, 10, 20, 40, 80, 125, and 160. We 
supplement this table by the table presented in Itoh et al. (1983) by using the same method 
for the case of r=l. In passing, it should be noted that those who wish to calculate the angle- 
dependent cross section should use the tables of the liquid structure factor iS'(A;) presented 
in Ichimaru, lyetomi, & Tanaka (1987) and also in Itoh et al. (1983). Constructing analytic 
fitting formulae for S{k) would be extremely difficult, since S{k) has oscillatory behavior. It 
would be best to use numerical tables rather than analytic fitting formulae for S{k). 

In Figure 2 we show the results of the calculation for < S* > in equation (15). A similar 
graph has been presented in Bruenn & Mezzacappa (1997). They have used Horowitz's 
(1997) fitting formula. 

In Figure 3 we compare our results with Horowitz's fitting formula. We find that 
Horowitz's results for e < 1.0 greatly differ from our results. It appears that this is due to 
the inaccurate treatment done by Horowitz to calcuate the small-A; behavior of the liquid 
structure factor S{k) that cannot be directly obtained from the Monte Carlo results because 
of the finite size of the simulation. We have used the correct small-A; (ajk < 1) behavior of 
the hquid structure factor (Ichimaru 1982; Itoh et al. 1983): 



S{k) = 



3r , 1 /(9P\ 



{ajky ksT \dn J , 



(21) 



1 f dP\ 

— = 0.73317 - 0.39890 F + 0.34141 V^'^ + 0.05484 V''^'^ . (22) 

ksT \dn J J, 

We find that the present more accurate calculation leads to a more dramatic reduction 
of the neutrino-nucleus scattering cross section for e < 1.0 than has been calculated by 
Horowitz (1997). Therefore, the ion-ion correlation effect on the neutrino- nucleus scattering 
in supernova cores will be more important than has hitherto been recognized. However, the 
consequence of the modification in the ion-ion correlation effect may not be so great as to 
drastically change the whole scenario of the supernova explosion, since the phase space for 
these low energy neutrinos is relatively small. 



- 7- 



4. ANALYTIC FITTING FORMULA 

In this section we present an accurate analytic fitting formula for < S* > (F, e) in order 
to facilitate the application of the numerical results obtained in the present paper. We have 
carried out the numerical calculations of < 5 > (F, e) for 1 < F < 160, 0.01 < e < 10.0. We 
express the analytic fitting formula for this range by 

< S >fit = — — — ' (23) 

X{j,0 = 1 - ^(e-l)'(7 + 2)'exp|-10(e-^-0.25)'-0.8(7-0.2)| 



-l(, + 2)-o.3exp{-6(e-^ + 0.2)^} 
+ ^(7 + 1-1)' exp |-130 + - 0.5)'| 



-^(7 + 2)'exp{-20(7 + 2)(e-0.25)^-0.8(7-0.2)} , (24) 

7 = (log.oF - 1.10206) , forl<F< 160 (25) 

^ = logioe , for 0.01 < e < 10.0. (26) 

The ratio of the present fitting formula value to our numerical result < S >fit / < S > is 
shown in Figure 4. It is readily seen that the accuracy of the fitting is better than 6% for 
all the parameter range considered. 



5. CONCLUDING REMARKS 

We have calculated the ion-ion correlation effect on the neutrino-nucleus scattering in 
supernova cores by using the liquid structure factor of the classical one-component plasma 
obtained by the improved hypernetted-chain scheme. We have found a dramatic reduction 
of the neutrino- nucleus scattering cross section due to the ion-ion correlation effect. We have 
compared the results of our calculation with those of Horowitz. We have found Horowitz's 
calculation greatly differs from ours for e = E^ai/hc < 1 because of his inaccurate treatment 
for the calculation of the liquid structure factor S{k) for small-A; values. Thererfore, we 
have found that the ion-ion correlation effect is more important than has hitherto been 
recognized. We have presented an accurate analytic fitting formula that summarizes our 
numerical results to facilitate the application. Supernova collapse simulations that take into 
account the present results are strongly awaited. 
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Fig. 1. — Physical conditions for the pure ^^Fe plasma. The present calculation is valid 
below the electron Fermi temperature Tp. Furthermore, the ions should be in the liquid 
state, r < 172. The lines of r=constant and r5=constant are also shown. 



10 



1.0 



10 



A ,„-2 
in 10 
V 



10 - 



10"" =. 



10"- 



-1.0 



-0.5 





























- ' ' ' ^^^^ .-^'V' 


r=i — 






r=2 - - - 






r=5 - - - 






r=io — 
r=20 — 
r=40 

r=80 — 
r=i25 — 
r=i60 



0.0 



log 8 = log 



EvOi 

he 



0.5 



1.0 



Fig. 2. — The angle- averaged ion-ion correlation factor < 5" > as a function of e = E^ai/tic 
for various values of F. 
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Fig. 3. — Comparison of the present calculation with Horowitz's fitting formula. 
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Fig. 4. — The ratio of the present fitting formula value to our numerical result < S >fit / < 
S>. 



